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Abstract 

The expression of survival factors for radiation damaged cells is em- 
pirical and based on probabilistic assumptions. We obtain it either from 
the maximum entropy principle for the classical Boltzmann-Gibbs entropy 
and/or from the Tsallis entropy. Empiric models are found to be particu- 
lar cases of the obtained expression. The survival factor exhibits a phase 
transition behaviour. This formulation supports different tissues grouped 
as universality classes. 
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1 Introduction 

Empiric laws are at the grounds of science. However, when a scientific discipline 
grows it is possible that the initially developed empiric laws become unable 
to explain all new experimental data obtained. Then engineering corrections 
and fitting coefficients emerge and usually explain those unexplained behaviors. 
When the experimental data grows enough even those corrections cannot explain 
the results; then a new set of corrections is looked for, found and a proper 
explanation for the new observed phenomena is provided. Eventually a large 
amount of anomalous behaviours is observed and the empiric law reaches its 
limit. It is impossible to go any further without getting deeper. 

At this point, a theoretical approach able to successfully explain the data 
is needed. Even if the new theory brings too complex expressions (like in the 
theory of relativity) , the limits of every approximation will be clear and it will 
be possible to reformulate the entire experimental approach. 

The tissue effect is an empirical concept widely used by oncologists to find, 
given a radiation dosage, the survival factor of tissue or tumor cells [H]. It is 
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widely used in clinical radiotherapy for both tumor and tissue cells and can be 
used to predict the outcome of a treatment. It has been developed for years and 
it is capable, as a concept, of gathering together several models of interaction 
between cells and ionizing radiation. Even when using the "single hit, single 
target" , the "single hit, multi-target", the "two hit, single target" or the linear- 
quadratic model [2| the powerful principle of grouping lethal events on a single 
adimensional quantity remains. 

The general expression used in radiobiology for fraction of survival cells is 
deduced assuming that lethal events per cell follow a Poisson distribution 
We think, besides, that the dependence between survival factor and tissue effect 
could be found from deeper laws based on first principles of nature. 

The maximum entropy principle is completely universal, well established and 
has an almost unlimited range of applications in physics, biology, demography, 
economy etc. In its modern formulation it says, "Given a model of probability 
distributions, choose the distribution with highest entropy." 0. It means that 
one should look for a distribution, consistent with the observed constraints, 
which maximizes the entropy. 



2 Radiobiology and extensivity 

According to the radiobiology linear model for the cell survival factor, the frac- 
tion of tissue killed by a radiation dose D is 

F d = l- exp[~aD] (1) 

Then F s = 1 — F4 is the cumulative probability of cell survival. From this 
follows that 

p (D) dD = aexp[-aD)dD (2) 

is the killed cell probability density (per dose unit). The cumulative probability 
fulfils the additive property: 

F s [D x + D 2 ] = F s [D x ] F s [D 2 ] (3) 

meaning that the effects of radiation are cumulative following an additive model. 
The additive variable here is the tissue effect, E — aD — — log [F s ] whose sum 
is 



E [D x +D 2 ]=E [Dx] + E [D 2 ] (4) 

However, the linear model is not accurate enough for a higher radiation 
dosage and the empiric experience shows that the tissue effect must be calculated 

as 

E = aD + f3D 2 (5) 

in what is called the LQ Model. The main problem involved here is that the 
probability function is not extensive. That is, the survival fraction for two doses 
is, 
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F s [£>i + D 2 ] = exp[-a(Di + D 2 ) - /3(D? + D\ + 2£> 1 £> 2 )] (6) 
and this implies, 

+ D 2 ] - F s [Di] • F s [£) 2 ] • exp^D^a] < F s [£>i] ■ ^s[A>] (7) 

This means that the survival fraction is lower if the radiation is applied in a 
single higher dose than if the same amount of radiation is applied in two doses. 
Or, in other words, a continuous radiation application kills more cells than a 
fractioned one. As E = — logFs we can arrive to 

E[D l + D 2 ] >E[D l ]+E[D 2 ] (8) 

meaning that E is a nonadditive variable. As a result of the nonlinear nature 
of E in this case, the superposition principle is not fulfilled. 

3 The classical approach 

Prior to reproduce the Tsallis entropy discussion for the nonextensive model, 
we will study the linear model extensive problem. This model fulfils that if the 
amount of absorbed radiation is unlimited, no single cell will survive, 

F s [£^oo]=0 (9) 

If the tissue effect is defined as proportional to the absorbed radiation E = 
olqD and p(E) is the probability density of killing a cell, then the fraction of 
killed cells for a tissue effect lower than E will be 

F d (E)= / p(x)dx (10) 
Jo 

and the survival fraction for every tissue effect lower than E will be 

/•OO 

F S (E) = / p(x)dx (11) 
Je 

Using p(E) it is possible to write a Boltzmann entropy functional, defined as 

s=J~dEm**^E) (12) 

Under the maximum entropy principle the expression of p(E) can be found 
if some assumptions are established. The completeness principle, 

/•OO 

/ p(E)dE = 1 (13) 
Jo 

and the mean value existence, 

/•OO 

/ p(E)EdE = (E) < oo (14) 
Jo 
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are demanded. 

Using the maximum entropy principle, the functional 

/■OO -1 nOO nOC 

/ p(E) log — — dE + a / p(E)dE + b / Eexp[-E]dE (15) 
Jo P{E) Jo Jo 

is built. Here a and b are the Lagrange multipliers and can be found following 
the imposed restrictions and maximizing the functional. It is straightforwardly 
obtained 

p{E) = ^y e_Tfy ( 16 ) 

the exponential distribution for the tissue effect. 
The survival probability of a single cell will be 

F s = e-Tk (17) 

We must note that ifTTj) is the experimentally proved and the normally used 
expression for the survival factor as a function of tissue effect and justified in the 
literature only through probabilistic arguments 0. We can take a — a>o/ (E) 
1 / (D) and the expression (fl7l) gets written in the known standard radiobiology 
form. 



4 The generalized approach 

Now that the extensive problem is solved, we will look for the non extensive 
solution in a similar way. To apply this principle we demand that there exists 
some amount of absorbed radiation A < oo, or its equivalent tissue effect, 
£1 = aoA, after which no cell survives, 

F s [Q}=0 (18) 

We will propose the use of Tsallis entropy 3, 



as the generalized entropy. We impose the conditions 



and 



(19) 



/ p(E)dE= / p{E)dE = l (20) 
Jo Jo 

/ p q (E)EdE = / p"(E)EdE = (E) < oo (21) 
Jo Jo 
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and apply the method of Lagrange multipliers. The functional 



l-J?p«(E)dE 
1-q 



p(E)dE + b q / dEp q (E)E 



(22) 



is maximized under those conditions and the values of 12, a q and b q are found, 



(E) q 
2-q 

(E), 



b q = - 



1-3 
1 



<l 



2-q 
(E) 



1- q 
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From where the probability density function 



(I 



1 - 



p(E) = 

is obtained. Then, the survival factor is 



1 



2 - 3 V <£) e 



(23) 
(24) 
(25) 

(26) 
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or using expression (J23j) , 
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(E) q 



E 



Fs(E) 
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(27) 



(28) 



for every E < fi. From 11231) it is not hard to see that when we deal with the 
extensive limit (q — > 1) then f2 — > oo as (£'} (J — > (E 1 ). 

Defining S = ao-D as in the previous subsection, we eventually get the 
expression for the survival factor of cells under radiation, 



Fs(D) 




VD < A 
VD > A 



(29) 



Notice that whereas (|28j) is expressed as function of the non measurable mag- 
nitudes ft and E, lf29|) expresses the survival factor in terms of the measurable 
quantities D and A. 

The value of A characterizes a critical point for cell survival probability 
and divides the phase plane in two sections with very different properties. For 
D < A probabilities of survival and death coexists but when D becomes equal to 
A, a phase transition takes place and no cell survives. This behavior resembles 
phase transition in ferromagnetics near the Curie point. 

Rescaling the radiation dose as D/A would allow to study the reaction of 
tissue cells under radiation in a more general way. 
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5 The q — > 1 limit 



The expression |27[) must include the particular cases corresponding to extensive 
systems. The cell survival probability limit (ft), fits into this principle. Expres- 
sion (f23]) gives the maximal tissue effect as a function of the g-mean value of E. 
(E) remains bounded for any value of q. When q tends to 1, (E) tends to (E) 
and Jl diverges in the extensive case. The obtained divergence is a trivial result 
that raises from Boltzmann formulation. However for q ^ 1 the divergence dis- 
appears in agreement with the experience. As shown in figure SI is big only 
if the system is extensive enough. If q is far from the unity, every cell will be 
annihilated with a finite amount of radiation. 

If the system is almost extensive (gal), and the radiation dosage is not too 
high, we can write the equation l|27p as the Taylor expression of the exponential 
function up to second order, 



F S (D) =exp 



(E) q J 2 2-q\ (E) q 

and comparing with ((5j) we find 



(30) 



° = *[WJ =A~ (31) 

(j - 2 2-q{(E)J ~2An-q (32) 

These expressions show that in the extensive case (q = 1), then a = 1/ (D) 
and (3 = recovering the linear model obtained from the Boltzmann formulation. 
This also allows to establish a correspondence between the known magnitudes 
a and (3 from the widely used LQ model and the newly defined parameters. 
Furthermore, the (3 parameter of the LQ model can be interpreted as a second 
order approximation of the non-extensive model. However, if D is high enough 
or A becomes a small quantity, the LQ model becomes useless to describe the 
survival factor. 



6 Comparison with experimental data 

To compare with experimental data we have selected some example plots in []]] 
and [H], where the survival factor is plotted as a function of D and separated 
curves are obtained for different radiation conditions. We used (|29|) to fit. As 
it is a function of D/A, all curves collapse to a single one. This means that A 
is the natural unit of D. Every plot of F s for a given tissue must converge to 
the same curve if D is expressed in the appropriate A units. 
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Figure [2] shows experimental data of a human melanoma under radiation for 
F s as a function of D/A at different dose rates. 

We can see that, at least for the represented ones, even when the dose rate 
affects the transition point A, the exponent j5f m p9| remains constant. 

This is also true for the kind of radiation. Figure [3] shows cell survival data, 
extracted from (jj, for stem cells under a beam of neutrons and electrons in 
different conditions. 

Though the cell survival factor (F s ) as a function of D depends on multiple 
factors, like the dose rate or the kind of radiation, the dependence of F s with 
D/A has a universal character. In this CclSG, clS the variable D is rescaled by A, 
the curves collapse and a universal behaviour emerges. The main factor then 
becomes the exponent, dependent on q, which in fact divides the cell behaviour 
under radiation in universality classes. The parallel with critical phenomena 
becomes apparent. 

7 Conclusions 

We have found a theoretical approach that puts the well known LQ model 
on physical grounds starting from first principles, rather than on probabilistic 
assumptions attempting to express the survival factors in terms of the tissue 
effect. Analyzing the used expressions we found that the tissue effect is not 
additive if defined through the LQ model. The non linearity of the tissue effect 
makes inapplicable the superposition principle, i.e., the tissue effect due to a 
continuous dose is higher than the corresponding to a fractioned one if the same 
amount of total energy is applied. 

The survival factor expression as a function of the absorbed energy per unit 
mass was derived from the maximum entropy principle following the Boltzmann 
entropy. Since the linear model is obtained from the Boltzmann entropy and 
it does not explain the known experimental data, so, a more general approach 
is introduced. Using the Tsallis g-entropy formulation a generalized expression 
for the survival factor is found. The extensive limit shows that the obtained 
expression is consistent with the known empiric laws. Besides, the empiric 
coefficients could be expressed in the new terminology and a new interpretation 
of its meaning is provided. 

The law obtained for the survival factor exhibits a phase transition behav- 
ior similar to second order ferromagnetic phase transitions where the imanated 
state corresponds to cell survival. A critical value of the absorbed energy marks 
the frontier between the non survivor cell region and the coexistence between 
survival and death probabilities. The transition between both regions occurs 
with a non integer critical exponent revealing a behavior similar to the ferro- 
magnetic phase transition at the Curie point. Rescaling every case with this 
coefficient allows to find the common shape for different experimental data be- 
longing to the same tissue. When compared with available experimental data, 
the phenomenon is shown as universal for a given tissue. Values of q would 
allow to group different tissues in universality classes. 
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Figure 1: Tissue effect limit (f2) as function of the extensivity parameter (q) for 
(E) q = 1. Even when in the plot (E) is assumed constant for every q, as the 
values of this parameter remains enclosed, the general qualitative behaviour, 
near q=l, must be close to the represented one. 
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Figure 2: Cell survival curves for a human melanoma irradiated at dose rates 
of IbOcGy/min (A = 28Gy), 7.6cGy/min (A = 42Gy) and l.QcGy/min (A = 
48Gy). The solid line is the survival factor as function of D/ A for q = 0.93. 
Data was extracted from Nl. 
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Figure 3: Survival curves for intestinal stem-cells. The radiations were neutrons 
(A = 39 Gy), electrons [high (A = 68 Gy) and low dose rate (A = 76 Gy)] and 
electrons under hypoxic conditions (A = 180 Gy). The solid line is the survival 
factor as a function of D/A for q = 0.97. Data was extracted from 
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